NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



f 



■ f* 


^K/tijA-CK-iu^i. oJ) ti-i-lCi Ck A iiNi'lt 
iUi4i/i Alibis hATii UN lilt n/tLi.AilVfc, liLAllrto Jr 
t LANh'I Ai nU:ii'iiLMu i;;»'I«Y PnCutJ 
(.1i^:_.oaLi JjiiV, -aUiici.) I, > j, /i 0‘4/Ilf AJ1 

<-bL L J .‘0 0 J/‘> 1 






0 lie i. .i i 
-71^vy 







lUASA 

Technical Memorandum 82094 


DIRECTIONALITY EFFECTS IN THE 
TRANSFER OF X-RAYS FROM 
A MAGNETIZED ATMOSPHERE: 
BEAM PULSE SHAPES 


P. M^szaros and S. Bonazzola 


(N AS A-ri*l-B2U^4) 
TtiS '1 HAN Si* SB Oi X 
Al’DlOSPiiSHE: HliAK 

lie ACJ/M'*' A01 


SiBECTlUNALITY Sii'Ec'iS IN 
-HA^S FHCM A MAGNET LZSU 
PULSE StiAtE (NASA) m 7 p 

CSCE GJS 


GJ/yo 


WU 1- JOOoJ 

LI nc L ct s 
3J2 7b 


FEBRUARY 1981 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 




Effect of a Finite Ionization Rate on the Radiative 
Heating of Outer Planet Atmospheric Entry Probes 


by 


H. F. Nelson 

Thermal Radiative Transfer Group 
Department of Mechanical and Aerospace Engineering 
University of Mlssourl-Rolla 
Rolla, Missouri 65401 
August 1981 


Prepared for 

National Aeronautics and Space Administration 
NASA Langley Research Center 
under Grant NASA NAGl-125 



CONTENTS 


Page 

ABSTRACT 4 

INTRODUCTION 5 

Present State of Knowledge 5 

SHOCKWAVE STRUCTURE 8 

NONEQUILIBRIUM MODELING 12 

Colllsional Pro- esses 12 

Radiation Processes 2 ^^ 

FINITE RATE IONIZATION — 18 

Leibowltz and Kuo Solution 20 

Tlwarl and Szema Solution 22 

Radiative Flux at the Body 26 

General Discussion 31 

SUMMARY AND CONCLUSIONS — 34 

REFERENCES — 36 

APPENDIX A: REVIEW OF NONEQUILIBRIUM SHOCKWAVE STRUCTURE A-1 

Relaxation Region Shock Structure In Atomic Gases A-1 

Ionization Relaxation A-3 

Nonequl librium Shockwave Structure In Hydrogen-Helium A-7 

Sliockwave Precursors — A-9 

Precursors In Atomic Gases A-9 

Precursors In Diatomic Gases A-12 

APPENDIX B: GENERAL FORMULATION OF RADIATION INTENSITY B-1 


1 


FIGURES 


Figure Title Page 

1 P “ 1 cm Hg, Trapped Radiation Mach 24 Shock Wave In Argon: 

Po • 2.285 X 10~’* gm/cm\ Tq * 16,828 K, Oq - 0.449. 9 

2 Schematic of Radiation Absorption Coefficient for a Two-Level 

atomic gas. 11 

3 Electron Temperature and Electron Concentration as a Function 

of Nondimenslonal Position across the Shock Layer, » 20 
represent^' Equilibrium flow, Q “ 0.95 represents Nonequlli- 
brluiii, 6(, » 1.7 cm. 21 

4 Temperature variation In Stagnation Shock Layer for Equili- 
brium and Finite Rate loniasatlon. 23 

5 Sp^icies Mass Fraction 'Mriation as a Function of Nondimen- 

slonal Position across the Stagnation Shock Layer. — — — 24 

6 Optical Thickness of the Stagnation Shock Layer as a Function 

of Radiation Energy. 25 

7 Equilibrium and Finite Rate Ionization Radiative Flux Toward 

the Body at the Stagnation Point. — 27 

A-1 Nonequilibrium Shock Wave Structure. A-2 

A-2 Atomic Transitions. A-7 

B-1 Radiation Geometry. B-11 


2 



TABLES 


Table Title PjOge 

1 Reaction Scheme and Rate Cons!;ants for Hydrogen-Helliim Non- 

equlllbrlum Conditions. 13 

2 Predictions of Equilibrium and Finite Rate. Radiative Fluxes. 19 

3 Population of the First Excited State of Atomic Hydrogen. 33 


3 



ABSTRACT 


The Influence of finite rate Ionization In the invlacld gas Just behind 
the stagnation shock wave on the radiative heating of ;;»robes entering the 
hydrogen-helium atmosphere of the major planets has been investigated. At the 
present time, there Is disagreement as to whether the radiative flux Increases 
or decreases relative to Its equilibrium value when finite rate Ionization Is 
considered. Lelbowltz and Kuo contend that the finite rate Ionization In the 
hydrogen gas just behind the shock wave reduces the radiative flux to the probe, 
whereas Tlwarl and Szema predict that It Increases the radiative flux. The 
radiation modeling used in the calculations of both pairs of these investigators 
has been reviewed. Tlwarl and Szema assumed a Boltzmann population at the local 
electron temperature for the electronic states of hydrogen. The electron tempera- 
ture just behind the shock is higher than the temperature for equilibrium flow 
field calculations. This increases the radiative source function and the ex- 
cited state populations, both of which Increase the radiative emission In the 
Baltner region of the spectrum. Because shock layers are optically thin in this 
spectral region, the increased emission directly Increases the radiative heat 
transfer to the probe. The finite rate ionization model developed by Lelbowltz 
involves the assumption that the excited state population is in equilibrium with 
the electrons at the electron temperature- It also Involves a finite excitation 
rate. This model predicts low excited state populations just behind the shock 
as compared to equilibrium chemistry solutions . The lower excited state popu- 
lations reduce the radiative emission just behind the shock and in turn decrease 
the radiative heating to values less than those predicted by equilibrium chemistry 
solutions. It Is concluded that finite rate ionization in the Inviscld region 
of the shock layer should reduce the cold wall radiative heating below the values 
predicted by equilibrium chemistry assumptions. 
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INTRODUCTION 


Radiation fortna a significant fraction of the heat transferred to a probe 

that enters the atmosphere of Jupiter or Satuvn. In reality, the design of the 

heat shield for such a probe is controlled by the radiation heat transfer. AK 

present, the United States r'oes not liave experimental facilities capable of 

producing the extreme thermodynamic conditions necessary to test the heat shield 

designs for Jupiter or Saturn ent>'y. Thus, the heat shields must be designed 

based upon analytical models of the flow field and heat transfer phenomena. 

The theoretical and numerical models must be as accurate as possible. In order 

(1 2 ) 

to Improve the heat transfer predictions, Lelbowltr and Kuo ’ and Tiwarl and 
(3 4 ) 

Szema ' Investigated the influence of finite rate ionization in the Invlscld 
part of the stagnation shock layer on radiative heating of probes during Jupiter 
and Saturn atmospheric entry. Both investigative teams considered cold wall in- 
viscid models for the probe shock layer. However, they reached opposite con- 
clusions as to how the ionization rate assumption affects the radiative transfer. 
(3 4 ) 

Tiwari and Szema ’ contended that the radiation heating at the probe stagna- 
tion point Increased, whereas Leibowltz and Kuo^^’"^ predicted that it decreased 
relative to the results obtained when the ionization rate was assumed to be in 
equilibrium at the local thermodynamic conditions. 

Present State of Knowledge 

The first studies of finite ionization rate effects on the radiative heat- 
ing of probes entering the atmospheres of Jupiter, Saturn, and Uranius were pre- 

( 1 , 2 ) ( 5 ) 

sented by leibowltz and Kuo and Howe . These studies Indicated that the 

assumption of a finite Ionization rate predicts significant reductions in the 

radiative heating as compared to equilibrium ionization predictions. Leibowltz 

( 2 ) 

and Kuo found that a nonequilibrium ionization layer with reduced radiative 
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emission exists just behind the stagnation point shock wave during part of the 

entry Into model atmospheres of Jupiter and Saturn. They predicted that the 

radiative heating rate, calculated using the finite Ionization rate assumption, 

could be up to A0% lower than if the heating rate were calculated by assuming 

that the plasma Ionization takes place In equilibrium. However, the reduction 

In total probe heating, which was calculated by Integrating both the radiative 

and convective heating rates over the atmospheric entry time, was less than 15%. 

The low reduction in total heat transfer occurred, because, when the finite 

ionization rate radiative fluxes were very low (compared to the equilibrium 

Ionization rate radiative fluKes), the heating -was largely due to convection. 

Leibowitz and Kuo^^^ calculated the radiative flux to the probe for both 

the equilibrium and finite Ionization rate assumptions by using the results of 

( 6 ) 

the hydrogen-helium radiation computations of Stlckford . They computed the 
Instantaneous values for the radiative flux at the thermodynamic conditions 
obtained from both the finite and equilibrium ionization rate assumptions at 
each point along the trajectory and Integrated the results to obtain the total 
radiative heating. 

(3 4) 

Tlwari and Szema ’ also considered the effect of finite Ionization rates 

on the radiative heating of probes entering the Jovian atmosphere. Their results 

showed that the assumption of a finite Ionization rate Increased the radiative 

heating rate by about 10% as compared to the heating rate obtained by assuming 

equilibrium Ionization. They used a fundamental fluid mechanics model and the 

(1 2 ) 

same ionization rate constants that Leibowitz and Kuo * used; however, they 
employed a detailed 58-step, spectral, absorption coefficient model for the hy- 
drogen-helium plasma to calculate the radiation transport. 

Zoby and Moss^^' conducted a preliminary Investigation of the thenrial en- 
vironment and heating rates for a probe entering the atmosphere of Saturn. They 
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used the numerical code developed by Tlwarl and Szema to determine the effect 
of a finite ionization rate In the Invlacid part of the stagnation shock layer 
on the heat transfer to the probe. Their results Indicate that the assumption 
of a finite Ionization rate would significantly increase the radiative heat 
transfer rate at the probe stagnation point as compared to equilibrium lonlza" 
tion and that the finite Ionization rate effects on radiative heating would be- 
come negligible away from the stagnation point. However, they did not allow the 
electrons and atoms to have different temperatures. 

These previous Invest igatlont, show that the assumption of a finite ioniza- 
tion rate influences the radiative transfer and, hence, the design of the heat 
shields of probes Intended for atmospheric entry of major planets. The non- 
equlllbrlum ionization laysr just behind the shock wave Influences the radiation 
transfer at or near the stagnation point. However, the radiative heating for 
the assumption of a finite ionization rate has been shown either to Increase 
or decrease relative to its value for equilibrium ionization for roughly the 
same physical conditions in two Independent investigations. This contradiction 
is investigated in the present study. 

The present study is limited to an Investigation of the underlying reasons 

(1 2 ) 

for the opposite conclusions reached in the work of Leibowltz and Kuo ’ and 

(3 4 ) 

Tlwarl and Szema ’ for nonviscous, hydrogen -helium shock waves with a cold 
nonblowing wall boundary condition at the probe heat shield. The effects of 
ablation and the ablation layer gases are not coi^sldered. The study is limited 
to the stagnation shock layer. Appendix A contains a survey of the present 
state of knowledge of nonequi librium shock wave structure. 
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SHOCK WAVE STRUCTURE 


This section contains a short Introduction to nonequilibrium shock wave 
structure. A detailed review Is Included In Appendix A. Figure 1 shows the 
shock wave structure for a typical atomic gas taken from reference 8. The 
figure Is for a Mach 24 shock wave In argon; however, It Is similar to that of 
a hydrogen-helluui shock wave, because the hydrogen dissociates almost Instantly 
behind the shock. The pressure In front of the argon shock Is 1 cm Hg, and the 
temperature Is 300 K. The parameters shown are all nondlmenslonal with tempera- 
tures, » T^/(16,800 K) , 1 - a, e, and density p ■ p/(. 000229 gra/cra^) . The 
reference value of the degree of ionization Is a_ “ 0.449. The parameter r 

^ O 

represents the ground state continuum optical thickness. The reference values 
of temperature, density, and are the equilibrium values behind the shock wave. 
Note that the nondlmenslonal parameters representing these variables all go to 
unity as the gas goes to equilibrium. The nondlmenslonal atom-atom and electron- 
atom Ionization rates are given by • 

The gas behind the shock wave Is In thermal nonequilibr ium, with the heavy 
particles much hotter than the electrons (f^ > f^) . As one moves into the re- 
laxation region behind the viscous shock, rapidly Increases to about its equi- 
librium value and then slowly Increases until it becomes equal to T and thermal 

a 

equilibrium occurs. 

Initially, just behind the shock wave, the degree of Ionization Is increased 
by atom-atom collisions. When the electron population Is sufficiently high, 
electron-atom collisions become the major producers of electrons. Each elec- 
tron produced via an electron-atom collision requires that the electron gas as 
a whole supply the necessary Ionization energy. Thus, the electron gas Is cooled 
by Ionizing collisions with the atom gas. At the same time, the electron gas re- 
ceives energy by elastic collisions with the atom gas and thereby reduces I • 

a 
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Figure 1. P = 1 cm Hg, Trapped Radiation, Mach 24 Shock Wave in Argon: 
= 2.285 X gm/ca^, Tn = 16,828 K, and a « 0.449. 



Throughout moat of the relaxation region, the energy gained by the electron 
gas through elastic collisions with the atom gas la almost balanced by the loss 
of energy through electron-atom Ionizing collisions. This keeps the electron 
temperature roughly constant throughout most of the relaxation region. 

The argon atom was modeled as If It had only two electronic utates: a 
ground state and an excited state as shown schematically In Figure 2. The 
single excited state of the mode , represents all the excited states of the atom 
Thus, the radiation model Involves only one bound-bound process and two bound- 
free processes. This type of twc-«cate model has been widely used In numerical 
shock wave structure studies. It generally predicts r-^tsults tliat agree well 
with experiment. 

Figure 1 shows that the ground-state continuum optical thickness, T , has 

8 

a value of about 8 at the end of the relaxation region. This implies that the 
ground-state continuum radiation emitted near the shock wave will not penetrate 
to the end of the relaxation region, because most of the radiation energy ab- 
sorbed within two optical path lengths of Its emission point. The excited 
state continuum optical thickness Is too small to be shown on Figure 1. Thus, 
excited state continuum radiation emitted just behind the shock can easily pene 
trate to the end of the relaxation region. However, the excited state source 
function Is very weak Just behind the shock, because there are not very many 
free electrons available to recombine to an excited atomic state, and In the 
process emit radiation. 
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KDNEQUILlBriUM MODELING 


Colllalon Frocesaea 

Leibowitz^^^ developed a reaction acheme deacrlblng the Important col- 
llslonal procesaea in hydrogen-helium Ionizing ahock waves that were mode3.ed 
after results obtained from studies of argon Ionization. Reactions are In- 
cluded for dissociation of molecular hydrogen, excitation of the electronic 
states of atomic hydrogen and helium, and Ionization of the hydrogen and 
helium by collisions with atoms and electrons. The reaction scheme and rate 
constants are listed In Table 1. 

In the reaction scheme, the atomic hydrogen and helium are assumed to have 
two bound states: a ground state and a firet excited state as shown schematically 
in Figure 2. The first excited state Is a representative state, which represents 
all the excited states. The two-state assumption is a good model when the ground 
and first excited states are separated by a large energy gap as compared to the 
energy difference between the first excited state and the Ionization energy. 

For these types of gases, the excited states and electrons have nearly the same 
energy. Thus, It is also assumed that the excited state and electron populations 
are in equilibrium with each other at the electron temperature. 

The cili:, 80 clation rate for molecular hydrogen was obtained from shock tube 
investigations . Hydrogen initially dissociates by collisions with Itself 
and helium and, subsequently, by collisions with atomic hydrogen, ions, and 
electrons. For the thermodynamic conditions of interest in outer planet atmos- 
pheric entry, dissociation is completed before ionization begins, so that the 
dissociation and ionization processes are uncoupled. 

The ionization of atomic hydrogen Is Initiated by atom-atom collisions, 
which produce atoms with electronically excited states. The excited atoms are 
ionized by additional collisions. The reaction scheme assumes that hydrogen 
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Table 1 Reaction scheme and rate constants for hydrogen- helium 
chemical noncquillbrlum conditions. (From Lelbowltz, 
Ref. 1) 


Reactions 


3 -1 -1 

Rate constants In cm sec mole 


1 . H + e H + 2e. 


- 2.27E13 Tg ' exp(-157,780/Tg) 


2 . He 4- e ^ He"^ + 2e 


k2 ^ 1.33E13 ' exp(285,160/Tg) 


3. H + eZH +e 


k 3 - A.11E13 exp(-116,010/Tg) 


A . He 4- e ^ He 4- e 


k^ “ 2.2AE13 Tg ^ exp(-232,030/Tg) 


5. H 4- H ? H 4 H 


kj - 6.20E10 T ' exp(-116,010/T) 


6. H 4- He 4 H 4- He 


k^ - 4.89E10 T exp(-116,010/T) 


7 f H 2 4- He H 4- H 4- He 


k^ “ A.33E18 [1 - exp(-15E8/T^)] 

exp(-52340/T) 


8 . H 2 + H 2 H 4- H 4 H, 


kg = 2.5 k^ 


9f H 2 + H ^ H 4 H 4 H 


kg ** lA.O k^ 


10. H, + H'*’ J H + H + 


ho - ''g 


11. H, +e;H+H+e 


'‘ii ■ S 


All temperatures are “K. 


Reaction rates k? and k^ modified from the expressions of Leibowitz by 
Tiwarl and Szema^^»^)_ 


n 




In the first excited state Is produced by collisions with hydrogen (reaction 5) 
and helium (reaction 6). Excitation to the first excited state Is the rate 
limiting step for Ionization. The excited atoms are rapidly Ionized by further 
collisions . 

As the result of atom-atom collisions, the population of free electrons 
Increases as one moves deeper into the shock layer. When the number of elec- 
trons reaches a certain magnitude, the electron-atom collisions become more 
Important than the atom-atom collisions. Electron-atom collielons produce 
hydrogen Ions by the same type of two-step process as the atom-atom collisions 
do. The electron is first excited by reactions 3 and 4 and then it Is Ionized. 

In summary, the reaction scheme developed by Lelbowltz Involves a two-step 
Ionization process In which the excitation to the first excited state Is the 
rate controlling reaction. In other words, the rate of excitation Is finite 
in the model and the first excited state Is not populated to Its Boltzmann 
value. Thus, the Lelbowltz model really allows for both nonequilibrium ex- 
cited state and electron populations. The finite rate Ionization model also 

includes the assumption that the excited state population Is in equilibrium 

(9) 

with the electron population at the electron temperature . Consequently, 
the model predicts low excited state populations just behind the shock wave, 
as compared to equilibrium chemistry solutions. 

Radiation Processes 

The general equations for the radiation absorption coefficient and source 
function are developed for a nonequilibrium atomic gas in Appendix B. The 
radiation absorption coefficient as given by Eq. 36, Appendix B, can be re- 
written for a two-level atom as follows: 

k^ “ “ Zj^^exp(-hv/kT^)l + n^(2)02£[^- “ Z2^exp(-hv/kTg) ] 

+ - Zj^2exp(-hv/kT^) (1) 
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The source function for a two-level atom from Eq. 39, Appendix B, becomes 


S - 


2hV 

2 


[ 


^”a^^^"^lf^lf na<2)a 2fZ2f)exp(-hv/kT^) n^(l)a^^Z^^exp(-hv/ kT^ 


-] 


( 2 ) 


where T. represents the electron temperature and T the atom temperature. 

The expression for and k^ as written above are for energies hv greater 

than the ground state Ionization edge. For energies (hv) less than the ground 

state Ionization edge, Is zero. Likewise, for energies (hv) less than the 

excited state Ionization edge, both and are zero. Also, a function 

of frequency, because it Includes the line shape. 

For the assumption that the excited states and electrons are in equilibrium 

at the electron temperature, one can write the expressions for ^2_2’ ^If* 

Z^, in terms of T and T . The expression for becomes 
2f e a ^ 12 


Z 


12 


n (2) 

_a 

n XD 

a' 



„^(1) gj 


exp(Xj^2^^’’a^’ 


(3) 


Wucre Xj^2 energy difference between the ground and excited state and 

where the bound electronic states are in equilibrium at the atom temperature 
when the atom is in equilibrium. Since the excited state is assumed to 

be in equilibrium with the electrons at T^, one can relate n^(2) to n^ by 


na<2) 


3/2 

(2nm kT ) g. 

2 -J- .xpC-X/Ha) . 


(A) 


where X 2 ionization energy for the excited state. If n^(l) is in equl 

librium with the electrons at the electron temperature, one has 
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n 


2(27ra kT g 

e e 


n*(l) 

a 


«1 


- exp(-Xi^/^T^) , 


(5) 


where ground stare ionization energy. Eliminating n^ from Eqs. (4) 

and (5), one can relate n^(2) to n^(l) and write 


n (2) n (1 )t g, 

rTTT TTrr^ ^ exp(-Xj^2^’^'^e^ * 


n (1) n (1) 
a a 


( 6 ) 


Thus, 2^2 written as 


^12 


n*U)T -X,,(T3 

-9 expl--12— 2- 


Te> 


n (1) 
a 


kT T 
a e 


1 . 


(7) 


where n (1) is given in terms of n and T by Eq. (5). Note that as the gas 
a .p e e 

approaches equilibrium Zj ^2 approaches unity. 

One can writ- as 


^ n*q) 

^If n^(l) n^d) 2 (Zron^kT^) 


h. exp(Xi/kTe ) 


( 8 ) 


where the free states and the bound states are in equilibrium at the electron 
temperature when the gas is in equilibrium. If one eliminates n^ by using Eq. 


(5), he obtains 


a T 
— fi- 


df n^(l) 
a 


(9) 


where n (1) is defined in terms of the local electron number density and the 
e 

electron temperature by Eq. (5). Note that approaches unity as the gas goes 
to equilibrium. 
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The parameter 1 b always unity, because the excited states and electrons 

are In equilibrium at the electron temperature. 

As a consequence of Eqs. (3) through (9), one can rewrite Eq. (1) for the 

absorption coefficient as 

n*(l) 

k » n (1)0,. {1 TiT^ exp(-hv/kT )} + n (2)o,,{l - exp(-hv/kT )} 

V a If n Ql) e a 2r e 

a 

"a^^^ -X (T - T ) 

+ ng(l)Oj^2^^ “ ®*P^kT^T — " — exp(-hv/kT^) } . (10) 

a' a e 


Equation (2) for the source function becomes 


. _ 2hV 

2 


3 


[n (1) O + n (2)0 1 exp(-hV/kT ) 

a T If ® 2f e 




+ n^(l) <J exp(-x (T^ - T^)/(kT^T^)l exp(-hv/kT^) 

2 ® 

\ 


( 11 ) 


As noted above, one must be careful in the application of these equations be- 
cause of the spectral behavior of 0^^^ and O^^* Note that as the gas approaches 
equilibrium, becomes the Planck function, and the absorption coefficient be- 
comes the product of the stimulated emission factor (1 - exp(-hv/kT)) and a 
term dependent upon the number densities of the bound states. The bound state 
population becomes Boltzmann distributed at equilibrium. 

Equations (10) and (11) show that the radiation source function and ab- 
sorption coefficient differ from their equilibrium values for the K2~He lonlza- 

(3 A) 

tlon model. Tlwari and Szema ’ used equilibrium values for the source func- 
tion and absorption coefficient. They evaluated k^ and at the electron 
temperature. As discussed below, this assumption for k^ and leads to in- 
creased radiative heating for the outer planet atmospheric entries as compared 
to a total equilibrium solution. 


17 



FINITE RATE IONIZATION 


The finite rate Ionization model for hydrogen-hellnm gas mixtures has been 
used In Investigations of the radiative heating of vehicles entering the atmo- 
sphere of Jupiter and Saturn. Table 2 lists a few of these Investigations and 
contains the ambient conditions, the probe nose radii, the shock wave stand-off 
distances, and the radiation fluxes at specific points In the entry trajectories 
for three cases. The solutions shown In Table 2 are for roughly the same entry 
conditions. Table 2 gives the comparison of the radiative flux for an equilib- 
rium solution and a finite Ionization rate solution attained from the same 

( 2 ) 

numerical code. The solutions obtained with the Lelbowltz code show that 
finite rate Ionization produces a decrease In the radiative heating rate, 
whereas those obtained with the Tlwarl and Szema^^’^’^^ code show that It 
produces an Increase in the radiative heating rate relative to equilibrium 
Ionization solutions. Note that In one of the cases, the equilibrium solution 

. fio) 

of Moss was used as a comparison case. 

The solutions of Moss^^^^ were obtained from an equilibrium viscous shock 
layer numerical code, which includes a detailed description of the radiative 
transport, equilibrium chemistry, and transport properties. The analysis al- 
lowed for ablation Injection at the heat shield and either laminar or turbulent 
flow. The set of equations describing the flow were obtained from the steady- 
state Navier-Stokes equations by retaining terms up to second order In the in- 
verse square root of the Reynolds number. This numerical code was used to pro- 
duce "benchmark" solutions for the design of the Galileo Probe heat shield. 

(3 4) 

Tlwarl and Szema ’ developed a viscous shock layer numerical code slml- 

( 10 ) 

lar to that of Moss to solve the steady-state shock layer equations. Their 
analysis Is restricted to nonblowing laminar flow fields about hyperboloid bodies. 
Equillbrltun and finite rate chemistry and thermal equilibrium or non-equilibrium 
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Table 2. Predictions of Equilibrium and Finite Rate Radiative Fluxes 



Conditions at Infinity 
in front of shock 

Probe 

Nose 

Radius 

Shock 

Stand-off 

distance 

radJ 

La t ion 
Elux 


“1 

km/sec 

T 

°K 

Poo - 
Kg/m^ 

% 

cm 

6 

cm 

mm 

• 

9non-eq 

Mw/m* 

Saturn Traj 
T “ 51.75 sec 
atm “ 89/11 

29.3 

83.9 

.000579 

31.1 

2.54 

40.0® 

71.7*’ 

(3) 

Saturn Entry 
Z » 115 Km 
atm * 85/15 

39.1 

145.0 

.000465 

23 

2.02 

422. 6*^ 

464.3^ 

(2) 

Saturn Entry 
T =• 18 sec 
atm * 73/27 

15 degree 

entry 



25 


37. 0^* 

20. 0^^ 


a. Moss, reference 10 | 

b. I’e “ Ta for non- equilibrium solution. Zoby and Moss, ref. 7 (using Tiwari and| 

Szema Code, ref. 3) J 

i 

c. Tiwari and Szema, ref. 4, pp. 161 | 

i 

d. Leibowitz and Kuo, ref. 2, Fig. 7 | 
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options are Included In the code. The equilibrium chemistry uses an approximate 
equation of state. The finite rate chemistry uses the rate constants and ioniza- 
tion model of Lelbowitz^^\ A 58-step spectral model for the absorption coef- 
ficient of the hydrogen -heliiim gas mixture was used to calculate the radiae^ion 
transport. 

( 2 ) 

Leibovitz and Kuo developed a simplified rapid numerical method to cal- 
culate a first order approximate solution for the flow field. They used cor- 
relations for the shock wave shape and stand-off distance and solved for the 
flow along stream tubes. Many details were omitted. The radiation calculations 
were uncoupled from the flow field calculations; however, they were correlated 
to the local thermodynamic state of the gas. A reaction similarity parameter, 
fi, which represents the ratio of the fluid mechanic residence time to reaction 
time, was used to correlate the electron distribution in the shock layer. When 
was small 1), the ionization rate was finite. When Q was large 20), 
the ionization was in equilibrium. 

Leibovitz and Kuo Solution 

( 2 ) 

Figure 3 shows a shock layer solution obtained by Leibovitz and Kuo for 
the electron temperature and the electron concentration for two values of at 
a shock layer position slightly displaced from the stagnation point. For the 
equilibrium solution, the electron temperature is approximately constant across 
the shock layer. In the finite rate case, T^ is high Just behind the shock and 
decreases to its equilibrium value near the body. For equilibrium flow, the 
electron concentration is roughly constant across the shock layer, and for finite 
rate ionization, it is very low near the shock and Increases to its equilibrium 
flow value near the body. In the Leibovitz finite rate ionization model, the 
excited state population is related to the electron concentration. In other 
words, it would be small just behind the shock and Increase to its equilibrium 
value near the body just as the electron concentration does. 
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Figure 3. Electron Temperature and Electron Concentration as a Function 
of Nondlmenslonal Position across the Shock Layer. U ■ 20 
represents Equilibrium Flow, ■ 0.95 represents Nonequilibrium 
Flow, 6 q ■ 1.7 cm. The Data Is for 0 ■ 0.20 Radians, where 
6 - 0 at the Stagnation Point. (From Ref. 2) 
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Tlwarl and Szema Solution 


(3 4) 

Tlviarl and Szeina ' used a detailed numerical analysis to calculate 
the flow field and to predict the radiation heating of outer planet atmo- 
spheric entry vehicles. They calculated the shock layer composition by using 
the finite rate ionization model developed by Leibowltz. Figures 4 and 5 show 
a typical solution for the electron and atom temperature and the composition 
at the stagnation point for both equilibrium ionization and finite ionization 
rate flow. The solution has the same general trends as that of Leibowltz and 
Kuo discussed above. 

Figure 4 shows that T and T for nonequilibrium flow are both greater than 

C ft 

the equilibrium temperature Just behind the shock, whereas, near the body the 
temperature is essentially the same for equilibrium and non-equlllbrlum flow. 
Thermal equilibrium occurs at about X) - 0.7 and the equilibrium and non-equl- 
llbrlum solutions yield about the same temperature for values of h less than 

0 . 2 . 

Figure 5 shows that the electron density Increases very slowly behind the 
shock for finite rate ionization (nonequilibrium) as compared to the equilibrium 
solution. This is an Important difference between the two cases when radiative 
emission is considered, because the electron population is very important in 
the radiation t -ansfer. Also, for finite rate ionization the atomic hydrogen 
population is much higher just behind the shock than the equilibrium hydrogen 
population . 

Figure 6 shows typical values of the shock layer optical thickness for the 
solutions shown in Table 2 as a function of photon energy hv. In general, the 
shock layer is optically thin in the spectral region from 1 to 9 ev. It is 
optically thick at the Lyman line at about 10 ev and for the Lyman continuum, 
at energies greater than 13.6 ev. The shock layer also becomes optically thick 
for free-free continuum radiation at very low values of photon energy. The 
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Figure 4. Temperature Variation in Stagnation Shock Layer 
for Equilibrium and Finite Rate Ionization. 
(From Ref. 4) 
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Figure 5. Species Mass Fraction Variation as a Function 

of Nondlmensional Position across the Stagnation 
Shock Layer. (From Ref. 4) 
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B’lgure 6. Optical Thickness of the Stagnation Shock Layer as a Function 
of Radiation Energy. 



lower part of Figure 6 shows the spectral position of the maxlaiuiD of the Planck 
function as a function of temperature. For the temperatures in the shock layer 
the maximum occurs In the spectral region where the hydrogen shock layer Is 
optically thin. Thus» when the radiation source function is the Planck function, 
most of the radiation emitted Just behind the shock can easily reach the body. 

Figure 7 shows the equilibrium and finite rate Ionization radiative heat 
flux toward the body as a function of nondlmenslonal position in the shock 
layer for the second case listed in Table 2. For the equilibrium Ionization 
case the majority of the radiative emission occurs relatively close to the 
body at about n “ 0.2. In the finite Ionization rate calculations the most 
intense radiative emission occurs near the shock wave at about n “ 0.9. This 
difference can be explained by considering the radiation source function and 
the optical thickness in the shock layer. 

Radiative Flux at the Body 

The radiative flux reaching the body involves an integral over the source 
function (see eqs. (40) and (41), Appendix B) . 

q^(0,V) = -2ir If S^(t^)exp[t^/y] dt^ dy. (12) 

-1 0 

Using the tangent-slab approximation, Eq. (12) can be integrated over y to yield 

q^(0,V) = -2n_/ S^(t^)E 2 (t^)dt^ (13) 

0 

where Exponential Integral of order 2. The source function 

S^(t^) and the absorption coefficient are functions of the electron tempera- 
ture and the electron and excited state populations as shown by Eqs. (10) and 
(11). When is small E 2 (x^ is approximately unity and the radiative flux can 
be written as 
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Figure 7. Equilibrium and Finite Rate Ionization 
Radiative Flux Toward the Body at the 
Stagnation Point. (From Ref. A) 
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qj^(O.V) 



s (t )dt . 

V V V 


This relation approximates the radiative flux at the body In the spectral region 
from 2 to 9 ev, where the shock layers are optically thin. Equation (14) shows 
that the radiative flux reaching the body Is the Integral of the source function 
over the optical radiating volume. It can be written In terms of the shock layer 
stand-off distance, 6, as 


qj^(O.V) 



S^(x)k^(X)dy. 


(15) 


Thus, the radiative flux at the body is directly related tc the radiation source 
function, the radiative absorption coefficient, and the shock wave stand-off 
distance when the shock layer is optically thin. 

Using Eq. (.11) one can rewrite Eq . (15) as 

^ 2hv 

2 
c 

0 


q^(0,v) * ^ f[n 


(1) a_ -t n (2)a„ Jexp(-hV/kT^) 
a ,j, If a 2f e 

e 


^ ^)«<P(-hv/kT^)] dy, 

a e 


(16) 


where it is assumed that the excited states and electrons are in equilibrium at 
Tg. Care must be taken in considering Eqs. (16) because 0^^^, and 
functions of the radiation energy. In the optically thin spectral region of the 
shock layer (approximately 2 to 9 ev) is zero. Also, one expects 0^2 6e 
small in this spectral region because it is mainly in the line wing. Therefore, 
one can write the radiative flux as 

3 

q^(0,v) = exp(-hv/kTg)dy. (17) 

•' c 
0 
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Thug, the radiative flux Is strongly dependent on and n^(2) . The value of 
0 ^^ for hydrogen from page 265 of Reference 11 is 

02J - 15.8 X 10"^®(hV2/hV)^ cm^ (18) 

where h \>2 is 3. A ev for hydrogen and hV is the photon energy in ev units. 

Tiwari and Szema assumed the source function to be the Planck function in 

their analysis. In addition, the absorption coefficient was determined using 

equilibrium populations for the atomic electronic states. In their computer 

code the absorption coefficient was calculated from the number densities for 

atomic hydrogen, ionized hydrogen, and the electrons at the local temperature. 

( 12 ) 

They used a 58-step spectral model developed by Sutton , which automatically 
assumed that the atomic electronic states were Boltzmann populated at the given 
temperature. For the thermal nonequilibrium, finite ionization rate cases they 
used the electron temperature to evaluate the source function and absorption 
coefficient. These assumptions are equivalent to writing the absorption coef- 
ficient from Eq . (10) as 


\ " exp(-hV/kTg) ] 


and the source function from Eq. (11) as 


3 -hv/kT^ 

S _ 2hy e e 

V " 2 -h\)/kT ’ 

^ 1 - e ^ 


(19) 


( 20 ) 


where Eqs. (19) and (20) only apply to the optically thin region of the shock 
layer. Further, Eq . (17) for the radiative flux at the body can be rewritten 
using Eqs. (19) and (20) as 


qR(0,V) 



0 


2hv^ 


c 


n^(2)o 


2t 


-hV/kT 

e , 

® dy . 


( 21 ) 
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Equation (21) points out the strong dependence of the radiative flux on the ex- 
cited state popvilatlon and the electron temperature. 

In the finite rate Ionization solutions of Tlwarl and Szema, Is high 
behind the shock as can be seen from Fig. 4. This Increases the exponential 
term of Eq. (21). High values of also Increase n^(2), because they assumed 
that the electronic state population Is Boltzmann distributed at T^. Thus, 
the finite rate Ionization solutions produce a large radiative flux. 

In the equilibrium case the temperature is smaller just behind the shock, 
as shown in Fig. 4. This reduces the value of the exponential term in Eq. (21). 
In addition, the number of hydrogen atoms Is much lower than it is for the finite 
ionization rate case (see Fig. 5). Both of these effects reduce the excited 
state population. Consequently, the contribution of the gas near the shock to 
the radiative flux is reduced, compared to the finite ionization case. Near the 
body Tg and n^(2) are about the same for equilibrium and finite rate ionization; 
therefore, the radiative flux from this part of the shock layer will be about 
the same for both assumptions. Thus, the overall magnitude of the equilibrium 
ionization radiative flux is less than the finite rate ionization radiative 
flux . 

The discussion to this point has been limited to the optically thin ap- 
proximation and the Balmer region of the spectrum. However, the conclusions 
are generally valid. As can be seen from Fig. 6 the ground state continuum 
radiation is optically thick so only gas near the body contributes to radiative 
heating of the body in that spectral region. The atomic line transitions are 
also optically thick near the line center, which means that the spectral regions 
near the line centers will contribute to radiative heating of the spacecraft 
only if the radiation is emitted close to the body. The free-free continuum 
transitions occur mainly in the spectral range below 2 ev. The free-free con- 
tinuum radiation can be important; however, at the temperatures of interest its 


contribution to the flux will be small, because of its small source function. 
Consequently, the radiative heating comes mainly from the Balmer region of the 
spectrum where the shock layer is optically thin and the above approximations 
hold . 

General Discussion 

(3 4) 

The physical problem Tiwarl and Szema ’ * solved is one in which the 
hydrogen atoms have a Boltzmann population distribution of their electronic 
states at the local electron temperature. The ionization rate was allowed to 
be finite by using the rate equations developed by Letbowitz^^\ They in- 
terpreted the excitation rates in the ionization model as given in Table 1, re- 
actions 3, 4, 5 and 6, as ionization rates. Also, they did not assume the ex- 
cited states and electrons to be in equilibrium at the electron temperature 
as Leibovitz did. Leibovitz used the excitation rate as the ionization rate 
limiting step in a tvo step ionization model. Once the electrons vere excited 
the number of electrons that became free was determined by the Saha equilibrium 
condition at the local electron temperature. 

In the Leibovitz ionization model the excited state is populated at a finite 
Thus, its population will lag its local equilibrium value. This implies 
that the excited state is not populated at its Boltzmann distribution value at 
the local temperature. Also, the excited state population and the electron 
population are assumed to be in equilibrium with each other at the electron 
temperature. Therefore, the excitation rate is not equal to the ionization 
rate. Both of these effects Influence the radiation transfer, because the ex- 
cited state population distribution influences both the absorption coefficient 
(Eq. 10)) and the radiation source function (Eq. (11)). In addition, Tiwari 
and Szema used the Planck function for the source function throughout the shock 
layer. As is pointed out above, the source function is equal to the Planck 
function only when the gas is in equilibrium. 
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Table 3 Illustrates the difference In the population of the first excited 
state of hydrogen for the assumption of a Boltzmann population distribution at 
the electron temperature and for the assumption of the excited state population 
being in equilibrium with the electron population at the electron temperature. 

The data for T^» and just behind the shock wave was obtained from the 
solution of Tlwarl and Szema shown in Figures 4 and 5. Note that for the Boltz- 
mann distribution the excited state population is much greater than it is for 
the assumption that the excited states and electrons exist in equilibrium to- 
gether at the local electron temperature. 

(3 4) 

The radiation absorption coefficient subroutine used by Tiwari and Szema ’ 
automatically used the Boltzmann distribution for the excited states population. 
This greatly increased the absorption coefficient compared to its value for the 
assumption that the excited state population is in equilibrium with the electron 
population at the electron temperature. This effect together with the high elec- 
tron temperature just behind the shock compared to the temperature for equilibrium 
solutions (resulting in a larger value for the radiation source function) in- 
creased the radiative emission and is the reason that Tiwari and Szema calculate 


increased radiation In the finite rate ionization case. 



TABLE 3. Population of the first excited state of Atomic Hydrogen assuming: 
1) Boltzmann distribution at Tg and 2) excited state in equilibrium 
with electrons at Tg. (Tlwri'i and Szema solution) 
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SUMMARY AND CONCLUSIONS 

Currently there Is disagreement as to whether the cold wall radiative 

heating for probes entering the atmospheres of the major planets Increases or 

decreases relative to its equilibrium value when finite rate ionization is con- 

( 2 ) 

oidered. Leibowitz and Kuo predicted that finite rate ionization in the hy- 
drogen gas Just behind the shock wave will reduce the radiative heating; how- 

(3 4) 

ever, Tiwarl and Szema ’ asserted that it will increase the radiation heating. 

The present study investigated the radiation modeling used in each of the 

above mentioned investigations. It is shown that the radiation transfer analy- 

(3 4 ) 

sis used by Tlwari and Szema ’ over predicts the radiation emission because 
they assumed a Boltzmann population distribution at the local electron tempera- 
ture (which is higher than the equilibrium solution temperature Just behind the 
shock) for the population of the electronic states of atomic hydrogen. This 
equilibrium population assumption changed the nonequilibrium radiative source 
function to the Planck function. In addition, it also influenced the av'^mic 
radiation absorption coefficient or both line and continuum radiation through 
the excited state population. Both of these effects increased the local radia- 
tion emission just behind the shock wave in the 2 to 9 ev region of the spectrum, 
where the shock layer is optically thin. This, in turn, increased the radiative 
heat transfer to the body as compared to equilibrium chemistry predictions. 

The more realistic and generally accepted assumption is that the excited 
states and electrons exist together in equilibrium at the local electron tempera- 
ture. This is the assumption used by Leibowltz^^^ in the development of the hy- 
drogen ionization model. In the finite rate ionization model, excitation to the 
first excited state is the rate limiting step. Once the electron is excited it 
is ionized according to equilibrium conditions at the local thermodynamic state. 
This Implies that the excited states are not populated to their boltzmann vsilues 
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at the local conditions. 


Tlwarl and Szema Interpreted the finite rate excitation reactions In the 
Ionization model as Ionization reactions. If they would have correctly used 
the finite rate Ionization model in their analysis they would have calculated 
smaller excited state populations and a smaller source function In the region 
just behind the shock wave. These effects would have reduced their finite rate 
Ionization radiative heating predictions to values below those predicted by as- 
suming equilibrium chemistry in the shock layer. 

In the near future the finite rate Ionization model as developed by Lei- 
bowitz should be correctly used In a detailed flow field and radiative heat 
transfer analysis like that of Tlwarl and Szema. The results of such an in- 
vestigation need to be available as "benchmark" solutions. Simplified solu- 
tions which make use of the optically thin shock layer assumption and are com- 
putationally fast need to be developed and compared to the "benchmark" solutions. 
Once they ire developed, the numerically fast, finite ionization rate schemes 
need to be upgraded to include the effects of ablation of the wall. The ablation 
products should he allowed to have finite chemistry and non-Boltzmann populations 
of their electronic and/or vibrational states. These effects should be included 
in predictions of the radiation heat transfer to probes entering the. atmospheres 
of the outer planets. The radiative heating predictions need to be as accurate 
as possible, because the heat shield mass on these missions are very critical. 

One does not want to over predict the size of the heat shield because of the 
weight penalty. Yet if one under predicts the size of the heat shield the probe 
may burn up before it completes its mission. 
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APPENDIX A: REVIEW OF NONEQUILIBRIUM SHOCK WA'^ STRUCTURE 


This section Is devoted to a survey of the conclusions and results of 

theoretical and experimental studies of the relaxation phenomena In shock waves 

and the shock wave precursors. Relaxation phenomena behind shock waves have 

(11 13-17) 

been discussed in reviews and monographs^ ’ , This section contains a 

discussion of the research conducted on the relaxation and precursors of non- 
equilibrium shock waves. This section does not contain an exhaustive account 
of all the research conducted on shock wave structure, but rather presents a 
discussion of the background experiments and theory leading to the present 
state of knowledge. 

Figure A-1 shows a schematic of the nonequilibrium shock wave structure 
for an atomic gas. It is typical of shock waves advancing at Mach numbers 
from 10 to 40, into a gas at 300 K. and 1 cm Hg pressure. Just behind the shock 
wave there is a collision dominated relaxation region, wherein the shocked gas 
relaxes to its thermal and chemical equilibrium state. In the equilibrium 
region the gas cools due to the emission of radiation. Some of the radiation 
emitted behind the shock wave propagates into the region in front of the shock 
wave, where it is absorbed by the cold gas to form the precursor. The degree 
of excitation and ionization in the precursor is coupled to the relaxation re- 
gion through the radiation transfer. 

Relaxation Region Shock Structure in Atomic Gases 
Strong sliock waves and high radiative heating rates that occur during the 
entry into the atmospheres of the outer planets have created a need to improve 
our understanding of the effects of the radiation, and collision processes in 
the relaxation process behind the shock wave. The hot gas radiative energy can 
propagate through the viscous shock wave and be partially absorbed by the cold 
gas in front of it. This absorbed radiation is referred to as trapped radiation. 
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Figure A-1. Nonequilibrium Shock Wave Structure. 



The absorption of this energy heats, excites and Ionizes the cold gas; there- 
fore, the shock wave propagates Into an excited and Ionized media. This pre- 
lonlzatlon and preexcitation can Influence the approach to equilibrium In the 
relaxation region behind the shock and the heat transfer to the entry vehicle. 
The radiation which Is not absorbed by the cold gas la lost, and Is referred 
to as radiation cooling. 

Heaslet and Baldwin^ have obtained solutions for shock waves propagating 

in a perfect grey gas in which all of the radiation is trapped. The nature of 

their solution depends on the amount of radiative heating compared to viscous 

dissipation. If radiative heating exceeds viscous heating, it Is possible to 

have a shock with no discontinuities. When radiative heating Is less tlian 

viscous heating, the shock has a discontinuity produced by viscous dissipation 

(19) 

which is embedded within a larger inviscid region. Cohen and Clarke and 
Chow^^^^ have shown tliat this is a valid model of the shock wave. If the pene- 
tration length of the radiation is large compared to the extent of the viscous 
region. Solutions for shock waves with combined radiation and viscosity with 

no restriction on the size of the penetration length have been obtained by 

( 21 ) ( 22 ) 

Traugott . Sen and Guess also solved the combined problem, but employed 
the Rosseland approximation to calculate the radiation transport in the shock 
wave . 


Radiation cooling effects have been observed in shock tube experiments by 
(23) (24) 

Petschek, et. al. , Redkoboradyl and Fedulov , and Oettlnger and Bershader 

(25) (26) (23) 

. McCheshney and Al-Attar^ and Petschek, et . al . , made calculations 


of radiation cooling assuming the shocked gas to be optically thin. Solutions 


including self-absorption effects have been obtained by Pomerantz 

( 28 ) 

Yoshlkawa and Chapman^ 


(27) 


and 


Ionization Relaxation 

Considerable information is available in the literature on the relaxation 
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of atomic gases behind strong shock waves. Most of the analysis has been done 
for shock tube ambient conditions of about 0.001 atmosphere pressure, 300 K and 


Mach numbers from 10 to 40. Experimental results have been obtained using argon 

rlcal 
(8,40-46) 


(25,29-35) (36,37) , , , (1,38,39) ^ , 

, xenon and hydrogen as the test gases. Numerical 


results for the ionization relaxation are available for shocks in argon 

(47,48) ^ , (49-53) ^ ^ (54,55) 

xenon , hydrogen and helium 

Ionization proceeds In a very complicated way behind shock waves In atomic 

(29) 

gases. Petschek and Byron showed that more than one ionization process was 
necessary to reach equilibrium ionization. If the electron concentration is 
sufficiently high, electron-atom collisions are the most probable ionization 
process. Petschek and Byron measured the electron-atom ionization rate. Their 
ionization rate was consistent with a two-step reaction scheme consisting of 
excitation followed by Ionization of the excited state. 


e + A A + e, 

A 

e + A *^A -fe-fe. 


Furthermore, their measurements showed that the electrons were not in thermal 
equilibrium with the atoms during the ionization process. This occurred be- 
cause the electrons gained energy by collisions with heavy particles less ef- 
ficiently than they lost energy by creating ions. Thus, ionization by electron- 
atom collisions resulted in a net loss of energy in the electron gas, which then 
became colder than the atom gas. Thus, the gas behind the shock wave can be 
assumed to exist as two separate gases: the electron gas and the heavy particle 
(atom) gas. The electrons establish equilibrium among themselves at the electron 
temperature while the atoms and ions establish equilibrium among themselves at 
the atom temperature. Consequently, thermal nonequilibrium will exist in the 
relaxation region. 

During the early stages of ionization the electron concentration is too 
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small to support Ionization by electron-atom collisions. Petschek and Byron 
speculated that reactions Involving atom-atom collisions, radiation, and con- 
taminants contributed to the Ionization In this region. Weymann^^^^ theorized 
that the most probable atom-atom reaction would be a two-step reaction similar 
to the two-step electron-atom reaction. 


A -f A A "♦•A, 

* -f 

A ■f'A"*'A + e + A. 


( 57 ) 

This reaction was experimentally verified by Harwell and Jahn . Impurity 
levels had to be reduced to a few parts per million before their effects be- 

( 58 ) 

came negligible. The experiment was further refined by Kelly who confirmed 

and improved tlie results of Harwell and Jahn. 

The effects of contaminant reactions have been investigated in xenon shock 

( 59 ) 

waves by Hacker and Bloomberg . They found a very complicated set of re- 
actions involving contaminant species and xenon atoms in various stages of ex- 
citation, and radiation. Contaminant reactions in argon were theoretically 
discussed by Morgan and Morrlson^^^^ . 

Radiation can be important in the ionization process too. Bibermann and 
Yakubov^^^^ have shown that line radiation from the region following the region 
of nonequilibrium ionization excites argon atoms in the nonequilibrium ioniza- 
tion region. Tlic excited atoms are eiaisily ionized and contribute additional 
electrons, which reduces the time to reach equilibrium. Radiation can affect 
the ionisation by creating electrons ahead of the shock front. Such precursor 

electrons have been observed in a re-entry by their effect on radar cross 

( 61 - 63 ) 


section 


and in shock tubes is described below. 


( 49 - 52 ) 

In a series of papers, Skalafurls investigated the structure of a 


shock wave in atomic hydrogen. His analysis was done for conditions that exist in 
stellar atmospheres ratlier than shock tube, or planetary atmospheric entry. He 


considered shock waves traveling at 30 to 70 kn/sec Into un-ionlzed atomic hy- 

-9 -10 

drogen at a temperature of 5000 K, for a range of densities from 10 to 10 

3 —5 —4 

gm/cro , which corresponds to a pressure range of 4 x 10 to 4 x 10 atm. 

The results show the same trends as those in argon, In that the Lyman radiation 
is Important in forming the precursor, the precursor raises the gas temperature 
behind the shock, the colllslonal relaxation region behind the shock Is optical- 
ly thin to Lyman and Balmer radiation and tliat the electron temperature Is less 
than the heavy particle temperature In the relaxation region. 

A general model for nonequilibrium shock wave structure In atomic gases, 
including both collisional and radiative processes, has evolved. It consists 
of a radiation Induced precursor, and embedded viscous shock, an inner col- 
lisional relaxation region and an outer radiation cooling region as shown in 
Figure A-1 . 

The relaxation mechanism to ionlzatlonal equilibrium occurs in two stages. 

The first is due to atom-atom collisions and the second Is due to electron-atom 
collisions. Thus, one expects a region just behind the shock where the electronic 
state excitation and the ionization are produced by atom-atom collisions. When 
the number of free electrons becomes significant the electron-atom collisions 
rapidly drive the excitation and ionization to equilibrium. Throughout most of 
the relaxation region the heavy particles are at a higher temperature tlmn the 
electrons because the electrons lose energy through ionization collisions faster 
tlmn they gain energy by elastic collisions with the atoms. The nonequilibrium 
gas usually reaches thermal equilibrium just prior to attaining ionizational 
equilibrium. 

A model has evolved for the excitation and subsequent ionization of atoms. 
Most atomic gases have a large energy difference between the ground electronic 
state and the first excited state compared to the energy difference between the 
first excited state and ionization as shown schematically in Fig. A-2 for an 
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Figure A-2 . Atomic Transitions. 
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atom with two electronic states. The electrons and excited electronic states are 


very close to each other in energy; consequently, they are assumed to exist 

together in equilibrium at the electron temperature. In an experimental study, 

( 9 ) 

Lelbowitz, et . al. found tliat most of the excited atomic electronic states 
were in equilibrium with the free electrons and were at the free electron tempera- 
ture. 

It has been shown experimentally that the rate of ionization is controlled 

by the rate of excitation between the ground state and the first excited state. 

Once the electron is excited to the first excited state it will rapidly ionize. 

Thus, the collisional ionization occurs as a two-step process as shown in Fig. 

A-2 . In atomic gas shock wave structure studies it is acceptable to assume 

that the excited state population is in equilibrium with the electron population 

at the electron temperature and that the ionization rate is controlled by the 

excitation rate between the ground and first excited state. 

Nonequilibrium Shock Wave Structure in Hydrogen-Helium 

( 381 

Belozerov and Measures^ ' theoretically and experimentally investigated the 
initial ionization process in strong shock waves in hydrogen. Their analysis 
indicated that the electron temperature and atom temperature were essentially 
equal throughout the relaxation region. This appears to be in error due to the 
use of an incorrect value of the elastic atom-electron collision cross sectlon^^^^ 
Lelbowitz^^^ obtained results that agree with Belozerov and Measures when he in- 
creased the elastic electron-atom cross section to the value used by them. Thus, 

T should be less than T through most of the relaxation region as occurs in 
e a 

argon shock wave structure. Nakagawa and Wisler^^^^ measured the electron tem- 
perature in hydrogen shock waves and compared excited electronic state populations 
with predictions. They found that the higher excited states tended to equilibrate 
with the electrons at the electron temperature much faster than the lower excited 
states . 
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( 19 ) 

Lelbowitz ’ studied ionization rates behind shock waves In hydrogen- 
helium mixtures by measuring hydrogen line and continuum radiation. He developed 
a reaction scheme which Included dissociation and a two-step excitatlon-lonlza- 
tlon meclianlsm for hydrogen Ionization by atom-atom and electron-atom collisions. 
He achieved good agreement between the numerical predictions and the experimental 
measurements. The excited states of hydrogen were assumed to be In equilibrium 
with the electrons at the electron temperature. The elextron temperature was 
found to be significantly lower than the atom temperature throughout most of 
the relaxation region. 


Shock Wave Precursors 

The presence of electrons In front of hypersonic shocks has been observed 

(66-80) , , . ^ u t (81-102) 

and theoretically investigated by several persons . These pre- 

cursors are important in atmospheric entry because they influence the propaga- 
tion of electromagnetic radiation in the vicinity of the entry body. This 
affects Identification of and/or communication with the vehicle. The precur- 
sor must be accounted for in interpreting the data a scientific probe gathers 

on atmospheric entry. Considerable effort has been spent determining the In- 

. . u (1-4,103) 

fluence of the precursor on heat transfer to atmospheric entry vehicles 

Precursor effects must be accounted for in the study of shock wave structure, 

since the precursor region is coupled to the relaxation region by the radiative 

transfer. The coupling of the precursor to the relaxation region has been stud- 


ied by several people 


(8,41,45,46,54,55,104,105) 


In addition, astrophysicists 


must understand precursors because they occur in the shock structure of variable 

(106) ^ (107) . ^ , (108) 

stars , interstellar gas , and certalyi solar events 

Precursors In Atomic Gases 

Throughout the literature controversy exists as to whether the precursor 
ionization observed in shock tubes is produced by electron diffusion from the 
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high pressure region behind the shock, or photolonlzatlon iron the ultroviolet 

radiation emitted in the region behind the shock. Photoemlsslon from the shock 

tube walls has also been considered as a possible source of the precursor. 

Many experiments have been undertaken to determine the mechanism for electron 

production. The earlier work leaned toward electron diffusion and photoemls- 
(66-69) 


slon 


; however, more recent experiments have conclusively shown that 


photolonlzatlon of the gas ahead of the shock is the main mechanism for pro- 

^ . (72-76) 

duclng precursors 

Precursor photolonlzatlon may occur either as a one-step process wherein 
the electron Is freed from the ground state of the cold gas In front of the 
shock by the absorption of ground state continuum radiation, or as a two-step 
process wherein the atom Is first excited by line radiation and then Ionized 
by excited state continuum radiation. These processes are shown schematically 
In Fig. A-2 . Lagar’Kov and Yakubov^^^\ Murty^^^\ Vulliet^^^\ Dobbins^^^^ 
and Nelson have investigated the Importance of line radiation in the pre- 
cursor and the coupling of line radiation and excited state continuum radiation 

in the precursor. Precursors caused by photolonlzatlon gaseous impurities have 

(74 

been investigated for trace amounts of hydrogen impurity In argon shock waves ’ 

90) (30) 

The electron temperature In the precursor has been measured 

The earliest precursor experiments were performed by Hollyer^^^^ in a metal 

pressure driven shock tube. He investigated Mach 9 shocks in argon at pressures 

less than one mm Hg, using Langmuir probes. Far ahead of the shock his probes 

collected only electrons. Thus, he concluded that photoemlsslon from the shock 

tube wall, due to radiation emitted from the shocked gas, was responsible for 

the precursor. 

Gloersen^^^^ experimented with xenon In a pressure driven pyrex shock tube 
at pressures from 0.75 to 4.0 mm Hg using probes wrapped around the outside of 
the tube. He observed two precursors, ofie traveling at about the speed of light 
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and one at the speed of the shock front. The fast precursor was attributed to 

photoemlssion from the shock tube walls, while the slow precursor was thought 

to be caused by colllslonal ionization, or photolonlzatlon of Impurities at the 

(67 7 7 78) 

shock front. Weymann, et. al. ’ ' used electrostatic and magnetic probes 

In a shock tube setup similar to that of Gloersen to Investigate argon shock 

waves at Mach numbers f,.'oro 8 to 12 and pressures from 2.5 to 10 mm Hg. They 

first concluded that the precursor was caused by diffusion of electrons from 

behind the shock wave because the electrostatic probe signals were negative 

and the magnetic probes indicated that electrons were moving ahead of the sliock. 

They reasoned that the large concentration gradients Just behind the shock 

and h:*i.gh thermal velocities of the electrons were responsible for the electron 

diffusion. However, a later set of experiments at lower impurity levels showed 

(79) 

that photoionization was the cause of the precursor. McRae and Leadon also 
found impurities to be the main source of precursor ionization in their experi- 
ments in xenon. 

(72) 

Zlvanovic investigated precursors with a ballistic range in which small 
copper spheres were fired up to 6,600 meters per second into air at pressures 
from 1 to 30 mm Hg. Electron density measurements were made utilizing a biased 
electrode plate, shielded by a screen, The bias established an electric field 
between the screen and the electrode plate. The field forced all the electrons 
which were photoemitted from the plate to return to the plate. Hence, only the 
electrons produced from diffusion and photolonizatiou contributed to the cur- 
rent through the plate, which was measured and related to the electron density. 
For a series of experiments a second grid was placed in front of the screen and 
biased so that it would repel all the diffusion electrons. Thus, these experi- 
ments measured only the electrons generated by photoionization between the grid 
and the screen. The measured current was the same in the experiments which mea- 
sured j'',ly photolonizetlon electrons and those which measured both diffusion and 
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photolonlzatlon electrons. Thus, It was concluded that the precursor was caused 
mainly by photolonlzatlon. 

I.ederman and Wllson^^^^ measured argon precursor electron densities with 
a microwave resonant cavity and an electrostatic probe. The Mach numbers varied 
between 10.5 and 13.5 and the argon pressures were between 0.76 and 15 mm Hg. 
Some of the experiments were performed with a "Venetian blind" In the shock tube. 
The 'Venetian blind" allowed electrons to diffuse through It and blocked the 
radiation. Their re-jults agreed with those of Zlvanovic In that photolonlza- 
tion, rather than electron diffusion, was the major cause of the precursor 
Ionization. 


Precursors in Diatomic Gases 

The literature Involving the Influence of precursors in diatomic gases is 

usually Involved In an engineering application. In most cases It Involves the 

effect of the precursor on blunt body heating for planetary entry conditions. 

(63,75) 

Precursors in diatomic gases have been observed and theoretically 

(93 9A 99 103) 

mode . '4 » » » ^ Much of the radiation energy emitted in the shock layer 

relaxation region is lost; however, when the photon energy is greater than the 
molecular dissociation energy, it is strongly absorbed in the precursor region. 
This absorbed energy dissociates and ionizes the gas. It also increases the 
gas temperature and pressure. The change in the free stream flow properties 
may change the flow characteristics of the shock layer and the heat transfer 
to the body . 

(75) 

Merniagen observed precursors in front of Lexan cylinders, launched from 
a light gas gun at velocities from 3900 to 5500 m/sec into air at a pressure of 
50 mn. Hg. The excited and ionized gases in shock layer around the model emitted 
infrared radiation which pre-excited the cold gas in front of the shock. The 


shock layer was populated by excited and ionized air, including contaminants, 
as well as ablation products from the Lexan. Lin'' presented evidence of a 
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’’photolonlxntlon lialo" surrounding the first manned Mercury orbital flight during 
Its re-entry. 

Radiative preheating of the cold gas by absorption of photons emitted from 
the hot shock layer was considered In most shock tube studies of shock wave 
structure by the late 1960's; however, it was usually neglected In the analysis 
of the blunt body shock layer problem. Yoshlkawa^ considered the precursor 
Influence on the blunt body shock layer for a gray gas and a linear relation- 
ship between the Planck emissive power and the gas enthalpy. Lasher and Wil- 
(103) 

son considered an uncoupled precursor, In which they used available blunt 

body shock layer solutions to obtain a first approximation of the radiative 

energy available to form the precursor. The precursor effects then preturbed 

the free stream cindltions and changed the shock layer solutions. Thus, their 

solutions were valid when the precursor effects were small. Lasher and Wilson 

Included the spectral variation of the radiation End radiation properties. 

(97) 

Liu investigated precursor effects on the heating of vehicles entering the 
earth's atmosphere in terms of the ratio of radiative to convective flux. The 
precursor influence was negligible ^hen the radiative to convective flux ratio 

was loss tlvm unity. 

(93) (94) 

Smith and Edwards calculated electron densities in the precursor 
of a shock preceding a blunt body for earth atmospheric entry. Smith used a 
simple microscopic model for air and included three-dimensional effects by 
modeling the body and shock layer as a point source. He calculated velocity, 
pressure, density and electron concentrations in the precursor region. Edwards 
used a more descriptive microscopic model for air and modeled the blunt body 
shock 'ayer as a constant tejnperature disk. A constant fraction of the radia- 
tion absorbed in the precursor was assumed to yield the electron concentration 
along the stagnation streamline. The influence of the precursor on the blunt 
body heating for lunar returns was predicted to be negligible. 
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In the past 10 years there has been considerable Interest in designing 

probes to enter the atmospheres of the outer planets, like Jupiter and Saturn. 

These atmospheres are thought to be composed of mainly hydrogen and helium. 

Thus, considerable effort has been put forth to Investigate the heating of 

vehicles entering hydrogen^hellum atmospheres. As a part of this effort Tlwarl 
(3-4) 

and Szema have Investigated the effect of precursors on the heating of 

vehicles traveling In hydrogen-helium atmospheres. They solved the radiation 
coupled precursor-shock layer problem taking Into account many details such 
as viscous flow, spectral details of the radiative transfer, thermal non- 
equlllbrlum In the shock layer, multi-dlmenslonal flow fields, finite chemical 
rates, etc. Their results generally show that the precursor effects should In- 
crease the heat transfer to the entry vehicle. 
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APPENDIX B: FORMULATION OF RADIATIVE INTENSITY 


This section of the Appendix contains a general formulation of radiation 
Intensity. The medium through vhlch the Incident radiation field passes Is as- 
sumed to be an atomic plasma. The plasma Is assumed to consist of ground and 
excited state atoms, Ions, and electrons. In foimulatlng the radiation Intensity, 
spontaneous and Induced emission and absorption will be considered. This 
model makes it possible to study the coupling of line radiation (bound-bound 
processes) with continuum radiation (bound-free processes) . 

The fundamental quantities of atomic radiation are the Einstein coef- 
ficients The Einstein coefficient, A^^, for spontaneous emission Is 

defined such that the probability per unit time for a spontaneous transition 
from upper state u to lower state S, with the emission of radiation in the 
solid angle dfi is 

A dfi, (1/sec) (1) 

ux. 

The total number of spontaneous transitions from state u to 5. per unit volume 
per unit time with the emission of radiation in dfl is 

n (u) A . dfi, (l/(cm^-sec)) (2) 

a ux, 

where n^(u) is the number density of the upper excited state. The energy 
emitted per transition is equal to hv^^ when the influence of the surrounding 
plasma is negligible; however, in real situations the transition can occur 
with the emission of almost any amount of energy. A certain value of emitted 
energy is determined by the emission line shape, (sec) which essentially 

represents the influence of the neighboring gas on the atom of interest. The 
total radiation intensity emitted in the frequency Interval V to V + dv per unit 
time is 

hv n (u) A „ d^dV (ergs/(cm^-sec)) . (3) 

a uJl V 
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Spontaneous emission also occurs when a free electron recombines to a 
bound state. In the model, there are several possible free-bound transitions: 
free to first excited (f to 2) and free to ground (f to 1), etc. 

The probability that an electron-ion encounter occurs with the free elec- 
tron losing energy, hv “ Xj^ emitting radiation within the solid angle 

is dil. After the encounter, the electron-ion pair becomes an atom 
with the electron In state 8-. For the model, Z can be any bound electronic 
level , 

Spontaneous emission cannot occur unless an electron and ion undergo 
an encounter. The number of collisions between electrons with energy in 
the range !i , t) + dri and ions per unit time per unit volume,, is 

A . n n f (n) dridi^ l/(cm^''sec.) (4) 

f £ el e 

where f^,(n) Is the energy distribution function for the electrons. For a 
Maxwellian distribution at the electron temperature T^ , 

f(il)dii = 2Ti(irkT^) ^'^^(n)^^^exp(-n/kTpdn . (5) 

The probability of a spontaneous emission, with the creation of an atom, 
is equal to tlie probability of an emission times the number of electron-ion 
collisions. The total radiation Intensity emitted In the frequency interval 
\’ to + d\' per unit volume per unit time in the solid angle diJ by all electron- 
ion spontaneous emissions Is 

hv A n n f (n) dndS^, (ergs/ (em'^-sec) ) (6) 

t < e i e 

where dri = hdv. This process emits continuum radiation; whereas, the pro- 
cess of Eq . (3) emits line radiation. 

The Einstein coefficient for absorption is defined such that in a 
radiation field of specific intensity, the probability per unit time for 

a transition from slate ? to u with the absorption of radiant energy in the 


B-2 


frequency range V to V + dv In the solid angle dfi Is 


B(Lu '^v (l/sec) (7) 

where is the absorption line shape. The total number of transitions from 
state £ to u per unit volume per unit time as a result of the absorption of 
radiation Is 

Bju (l/cm^-sec)) (8) 

The total amount of radiation absorbed from the radiation field In the fre- 
quency Interval V to v + dV per unit volume per unit time in the solid angle 

dP. is 

hv n^(0 dftdv, (ergs/ (cm^-sec) ) (9) 

This process represents a bound-bound process. The line shape, i*'’ used 
because the probability of the absorption of energy is related to the magnitude 
of the energy, through the line sliape. 

The bound-free process must also be considered. Equation (9) can be used 
to represent the bound-free process if the subscript u is changed to f. The 
continuum absorption can then be written as 

hv Bj^^ dfidv. (ergs/(cm^-sec)) (10) 

where k' represents any bound electronic state. 

Emission of radiation may also be -induced by the radiation field. The 
Einstein coefficient for Induced emission is defined such that in a radiation 
field of specific Intensity, the probability per unit time for a transition 

from state u to ? as a result of Induced emission in the frequency range V to 
+ dv and in the solid angle dSJ is 

B „ I dildv. (11) 

u!l V 
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The total number of transitions from state u to S, per unit volume per unit 
time as a result of Induced emission Is 

dfi, (l/cm^-sec)) (12) 

The total amount of radiation emitted In the frequency Interval V to V + dv 
In the solid angle dii per unit volume per unit time by Induced emission Is 
then 

hv n^(u) dfidv, (ergs/ (cm^-sec) ) (13) 

where again the line shape, Is introduced because It is the probability 
density function of the magnitude of the emitted energy. 

Equation (13) represents bound-bound induced emission. The Induced emis- 
sion which results in continuum radiation must also be considered. Induced 
emission, when the initial states are free states, depends upon the number of 
collisions between the electrons and ions. Using Eq . (4) and Eq . (13), one can 
write an expression for the continuum induced emission. The radiation emitted 
In the frequency Interval, v to V + dv, in the solid angle dfi per unit volume 
per unit time by free-bound induced emission is 

hv dndil, (ergs/ (cm^-sec) ) (14) 

where dii = hd\'. During the process, the electron loses energy, hv ■= Xjj + U, 
while tile ion gains energy, the recombination is at the dth level of the 

atom, as the electron-ion pair create an atom in state i. . 

The rate of cliange of the radiation intensity at a point within the solid 
angle, dSl, in the frequency range, v to V + dv, in the direction s is the sum 
of the contributions (both negative and positive) given by Eqs. (3), (6), (9), 
(10), (13), am! (14): 
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( 15 ) 


dvdi) - {-(n„a) B^^ f^(n) 

- %<“> ‘ui ♦v' \l K 

e 1 ti!, e 


The processes represented on the right-hand of Eq. (15) In the term multiplied 
by are respectively continuum absorption, Eq. (10), line absorption, Eq. 

(9), continuum induced emission, Eq . (lA) , and line Induced emission, Eq . (13), 
The second term, which Is multiplied by hV represents line spontaneous emission, 
from Eq . (3) and continuum spontaneous emission, from Eq . (6), respectively. 

Equation (15) represents the general form of the radiation transfer equa- 
tion. It can be rewritten as 


dl 


\) 


n (u) B n (J) 

+ n (V B (J) [1 - hvl 

a *t.u V n„(£) B„ <}) V 

a Jiu V 


(16) 
A 


n^(u)A „ i|; n n, .. o 

-f- n (£) B 4) — ~ hV -f n B f (n) hV , 

a Jtu V n^(5.)B^^ (|.^ ?, U n^W e 


il = 1,2, . . . L-1 
u = 2 , 3 , . . . L 

where L is the total number of bound electronic states. In what follows the 
absorption and emission line slopes will be taken as the «ame, 

The relationships between the Einstein coefficients can be found by con- 
sidering the equilibrium case. If the gas is in equllibriiim, = 0 and I^^ 

is given by the Planck function 




2hv~^ 1 

2 >hVv 1 

c exp(-j^; - 1 


( 17 ) 


Considering only the line transitions, those terms involving in Eq . (16), 
one can write at equilibrium 
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1 


( 18 ) 


a iu 


a jHu 


1 - 


However, at equilibrium 


n*(u) 
_ili 

n*(jt) 

a 


g„ 

— exp(- X „/kT ) 

go « 


(19) 


where and g are the degeneracies of the states H and u respectively, 
is the atom temperature, and is the energy difference between the states. 
To force Eq . (18) to agree with the Planck function one must have 


'^ut H 2hv^ 


and 


In 

®nu ' ®u 


(20a) 


or that 




r>] 

Ln*(u)l 


hv 


exp(- 


EdL'i 


kT 


un 

*nu 


Ln*(u). 


exp( •-^) 

kT^ 


(20b) 


Then one can write I., as 


2hv 


\) 

3 [n (u)/n (n)] [n (2.)/n (u)]^i exp(- hv„ /kT ) 

^ ^ X/U 3 


c. 1 - [n (u)/n (g,)] [n (g,)/n (u) ]- exp(- hv /kT ) 
a a a a I _ x,u a 


(20c) 


which approaches the Planck function as the gas goes to equilibrium. This 
then defines the three Einstein coefficients, if one of them is known. 

When the continuum portion of the intensity is considered, one must define 
the Einstein coefficients as 




( 21 ) 


where is the degeneracy of the free election states with energies between 


n and n + dn; 


(111) 


2(4tt) (mg) 


3/2 


ry ^ 1/2 ^ 

(2n) dn. 


( 22 ) 
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The factor Is the degeneracy of the Hth bound state, and is the degen- 
eracy of the ion. The Ion Is normally in Its ground state. In Eq. (22), 

3 

has units of 1/cm , so that it represents the degeneracy per unit volume. One 


n„n, B 

must consider the relation — j-pr f_(n) "T 


fi 


where 


Af 


f^(n) dn 


2n 


(TrkTp) 


3/2 


(n)^^^exp(- n/kT ) dn 


(23) 


represents the fraction of electrons with energies between n and H + dn, where 


ri 




•f hV. The relation can be written as 

n/kT^) 

or : :^™T^ ^3/2 


2g^ (ZrmigkT^) 


using Eqs. (21), (22), and (23). At equilibrium the ratio of the population 

(4) 

of the electrons and ions to that of bound state i is given by Saha equation 

3/2 

Zi. exp(-X£/kT^) ^ 24 ) 


•k A 

n n 

n (0 T 
e e 


= 2 


Thus, the relation can be written as 


n n 
e i 

n ( ?-) 
a 


n*(£) 

[-j— jt"] exp(-hv/kT^) 

n , n T 
lee 


(25) 


where the definition of n was used (n = Xjj^ + hV) . Using Eq. (25), one can 

write the continuum intensity as 

n n n (JJ.) 
e i r a 


I = 
\i 




n n n*(£) 

- • ^7 " oV exp(-hv/kT ) 

"a^^^ "e“i T 

e 


(26) 


Now as the gas goes to equilibrium and using equation (21) for 

I B , the Planck function. Note that Eq. (26) is a function of the electron 

\' V 
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temperature T^. 

Using the definitions given by Eqs. (20), (21), and (25), one can rewrite 
Eq, (16) as 

dl 

- Z^^exp(-hv/kl_^)] + 

• (1 - Z oxp(-hv/kT)J)hvI^ (27) 

3 

+ {n^(£)Bj,^Z^^exp(-hv/kT^)hv + n^(Jl)B^^<|)^Z^^exp(-hWkT)hv}-^ 


where 


and 




n n 
e 1 

n^(« 


n <u) 
a 

n (i) 
a 



( 28 ) 


(29) 


and ? Is summed over all bound levels £ = 1,2. . . L and u is summed over 


levels -f 1, f 4- 2 , . . . L. If the gas is in equilibrium, and Z^^ go to 

unity and dl /ds goes to zero. The general form of the equation of radiative 


transfer Is 


V 

( 110 ) 


dl^ 

ds 


- - Vv + <=v 


(30) 


where is the absorption coefficient, and is the emission function. From 
Eq . (27), one sees that 


K = ■ Zj_fexp(-hv/kT^)] 

+ n (^.)B <|) [1 - Z exp(-hv/kT) ] )hV 

a Jcu V £u 

and 

‘'v ' l"3 0t)B„Zifexp(-hv/kT^) + n^(S,)B^_^*^Z,^exp(-hv/kT)h\. 

c 


(31) 

(32) 
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In the limit of equilibrium, and for the lines agree with relations given 

by Jefferies^^^^^ when the absorption and emission profiles are the same. The 
radiation transfer equation Is simplified by converting the Independent vari- 


able from s to T , the optical depth. 




( 33 ) 


The equation of transfer can then be written as 

"v 

T?; - 'v - s. 


( 34 ) 


where ^ known as the source function, and y » ps, where y Is 


the cosine of the angle between the general direction of radiation Intensity 
propagation s and the coordinate y . The source function becomes 


hv 


hv 


2hv' 


hv 


hv 


( 35 ) 




The radiative transfer can now be evaluated, if the line shape, c})^» is known 
and the Einstein coefficients, and can be determined. B^^ and B^^^ 

are functions of the particular gas of interest. Note that in the limit of 

equilibrium, reduces to the Planck function. Note also that the units of 

2 2 
are cm /erg and that the units of are cm /(erg-sec) because must 

2 

have units of cm /erg. 

The absorption coefficient can be written in terms of the cross sections 


as 




jlV 

kf 


U 


_ JlV 

e) + n (^)cr„ (1 - Z e" kT) 
a )6u £u 


( 36 ) 


where the cross sections are defined as 




( 37 ) 
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and 




( 38 ) 


In terms of the cross sections, the source function becomes 


2 hv _ hv 

„„(J.)a„(l - z„o' “e, + n^dOaJl - “) 


(39) 


The cross sections are functions of the gas and the various transitions. The 
general solution of the equation of transfer for the geometry shown in Fig. B-1 


is 


( 110 ) 


T (t,,M) “ 

V V 


I 1 r" 


<vV'" ‘“v 


- (u<0) 


J V IJ 


-T /M 

(0,u)e ^ + 

V 


(u>0) 


(40) 


where (tJ ,li) is the intensity incident on the boundary at , and 

(0,)j) Js the Incident intensity on the boundary at zero. The contribution 
to the radiative flux in the frequency range V to V + dv is 


-1 


qj^(x^^,V) » 2TT / Ul^(T^,,.d) dM. 


V V 


(41) 
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Figure B-1 . Radiation Geometry. 
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